Application of screened hybrid functionals to the bulk transition metals Rh, Pd, and 
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We present the results of calculations on bulk transition metals Rh, Pd, and Pt using the screened 
hybrid functional YS-PBEO [F. Tran and P. Blaha, Phys. Rev. B 83, 235118 (2011)]. The results for 
the equilibrium geometry are compared with those obtained from (semi)local functionals, namely, 
the local density approximation and the generalized gradient approximation PBE of Perdew et al. 
[J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)]. It is shown that 
the screened hybrid functional yields more accurate equilibrium geometry than PBE, but, overall, 
it is not more accurate than LDA. However, in contradiction with experiment, we find that the 
screened hybrid functional favors a ferromagnetic state as the ground state for all three transition 
metals. Therefore, the use of hybrid functionals for, e.g., the study of catalytically active systems 
with correlated oxides on a metal support is questionable. 

PACS numbers: 71.15.Ap, 71. 15. Mb, 71.15.Nc 
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I. INTRODUCTION 

The Kohn-Sham (KS) version of density functional 
theory (DFT) is the most used quantum method for the 
calculations of properties of matter^ however, since 
we do not know the exact exchange-correlation func- 
tional, one needs to make a choice for an approxi- 
mate form for practical calculations (see Ref. y for 
a recent review). Many different types of approxi- 
mate exchange-correlation functionals have been pro- 
posed. The two most commonly used approximations 
in solid-state physics are the local density approximation 
(LDA) 2 . and generalized gradient approximation (GGA), 
which lead, in most cases, to reasonably accurate results 
for the structural parameters in particular (see, e.g., Refs. 
iandi). 

However, a well-known problem of LDA and GGA, 
is their inability to yield accurate excited-states prop- 
erties, but actually this is a more general problem which 
has its roots in the KS method itself when an orbital- 
independent potential is used (see Ref. @). The most 
common way to get more accurate excited-state prop- 
erties is to work in the framework of the generalized 
KS scheme^ by using functionals which lead to orbital- 
dependent potentials like DFT+J7 8 or the hybrids. 9 
DFT+U is computationally as cheap as a LDA/GGA, 
but can be applied only to well localized strongly cor- 
related electrons (typically 3d or 4/ electrons). The 
hybrid functionals, which mix semilocal (i.e., LDA or 
GGA) and Hartree-Fock (HF) exchange, can be used 
more widely, however for solids, they lead to calculations 
which are one or two orders of magnitude more expensive 
than LDA/GGA and their applications to metals can be 
problematic f^ - 

A way to reduce partially the problems encountered 
with hybrid functionals for solids is to get rid of the 
long-range HF interaction by screening the Coulomb po- 
tential with, e.g., an exponentiate or error function^ 



Actually, screened hybrid functionals are becoming more 
and more popular and have shown to be successful for 
structural and electronic properties on a large variety of 
solids (see Refs. fl3H2(il for extensive tests) including diffi- 
cult cases like transition-metal (TM) oxides and nitrides, 
where band gaps and magnetism are strongly enlarged in 
agreement with experiment (see, e.g., Refs. I21H27T ). In 
the case of a GGA-based screened hybrid functional, the 
functional reads 



E^ c — E. 
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where i?| R " HF and E^ K ~ GGA are the short-range (SR) 
parts of the HF and GGA exchange functionals, respec- 
tively, and a x (€ [0, 1]) is the fraction of exact exchange. 
The most known screened hybrid functional is the one 
proposed by Heyd et al. (HSE), 12 which is based on the 
GGA functional PBE of Perdew et al& (like the un- 
screened hybrid functional PBEO— >22) and uses the error 
function for the splitting of the Coulomb potential. Most 
of the time, HSE is used with a x = 0.25. Another known 
screened hybrid functional is the early one proposed by 
Bylander and Kleinman, 11 which is based on the LDA 
functional and uses a screening parameter that is cal- 
culated from the average value of the valence electron 
density (in HSE the screening parameter is fixed to a 
"universal" value). 

The trends of the HSE06 31 functional is to be very 
accurate for band gaps smaller than ~ 5 eV, but to un- 
derestimate larger band gaps. Actually, the larger the 
band gap, the more the underestimation (in percentage) 
becomes severe (see Refs. [H and [32]). Usually, GGA 
functionals underestimate the atomic magnetic moments 
in strongly correlated systems, while the values with 
HSE06 (or any other hybrid functional) are increased and 
hence in better agreement with experiment. 24 Concern- 
ing the lattice constants ; 19 ' 20 HSE06 is more accurate 
than PBE, but not as accurate as the GGA functionals 
which were especially designed for solids like WC— or 



2 



PBEsol. 34 However, we note that the values of the lat- 
tice constants depend strongly on the GGA functional 
on which the hybrid functional is based, while the band 
gaps are rather insensitive (see Ref. Il9). Actually, re- 
cent studies reported calculations on solids using hybrid 
functionals based either on the WG^ or PBEsol GGA 
functional, and in both cases it was shown that the hy- 
brid functional benefits from the good performances of 
the underlying GGA (WC or PBEsol) for geometrical 
parameters. We also mention that in Ref. [3(| a middle- 
range hybrid functional was shown to perform very well 
on semiconductors for the lattice constant and band gap. 
Recently (Ref. |37l ), a screened hybrid functional based 
on the GGA PBE functional, but using the exponen- 
tial function for the splitting of the Coulomb potential 
(known as the Yukawa potential), was implemented into 
the WIEN2k code.— The choice of using the exponential 
function instead of the error function was made for tech- 
nical reasons (e.g., integrals calculated more easily), but 
is of very little importance for applications. Indeed, it 
was shown that using the screening parameter A = 0.165 
bohr -1 (A was kept fixed at this value for the present 

work) in the Yukawa potential e~ A l r ~ r 1/ |r — r'| leads 
to results which are close to the results obtained with 
HSE06 31 for which ^ = 0.11 bohr -1 in the SR potential 
erfc (/i |r — r'|) / |r — r'| . The functional was called YS- 
PBEO (where YS stands for Yukawa screened). 

In this work, we present the results of calculations 
obtained with the screened hybrid functional YS-PBEO 
on the face-centered cubic 4d-TM Rh and Pd and 5d- 
TM Pt. Usually, the semilocal functionals describe fairly 
well the electronic and magnetic properties of the (non- 
magnetic 3d, 4d, and 5c? TM, although in particular for 
GGA functionals a tendency to overestimate magnetism 
can been founder— This is in contrast to the TM ox- 
ides (in particular the 3d like MnO, NiO,— or the un- 
doped cuprates^), where magnetism is underestimated 
and semilocal functionals give qualitatively wrong results 
due to the localized nature of the d-electrons, while in 
the pure TM the electrons are itinerant (i.e., more free- 
electron like). 



II. RESULTS 

The calculations were done with the WIEN2k code^ 
which is based on the full-potential (linearized) aug- 
mented plane-wave and local orbitals method to solve 
the KS equations^ For comparison purposes, calcula- 
tions with the semilocal functionals LDA (using PW92 47 
for correlation) and the GGA PBE2& were also done. The 
semilocal calculations were done with a 24 x 24 x 24 k- 
mesh for the integrations of the Brillouin zone, while a 
12x12x12 k-mesh was used for the hybrid functional YS- 
PBEO to make the computational time reasonable. We 
checked that using such a k-mesh leads to results which 
are accurate enough for our purposes. The size of the 
basis set is determined by the product R M fK max of the 



TABLE I: Equilibrium lattice constant ao (in A), bulk 
modulus Bo (in GPa), and the total-energy difference (in 
mRy/atom) between the minima of the FM and NM phases. 
The experimental values are for T — K and are corrected 
for the zero-point anharmonic expansion.— 



Method a$ M 



FM 

a 



nNM oFM 
-DO &0 



Rh 

LDA 3.757 316 
PBE 3.831 257 
YS-PBEO 3.799 3.874 279 
Expt. 3.786 277 



226 



0.59 



Pd 

LDA 3.840 228 
PBE 3.943 3.946 168 
YS-PBEO 3.935 3.938 165 
Expt. 3.876 187 



165 
168 



0.15 
3.96 



Pt 

LDA 3.896 305 
PBE 3.971 246 
YS-PBEO 3.948 3.959 265 
Expt. 3.917 286 



261 



3.51 



smallest of the atomic sphere radii (-Rmt) and the plane 
wave cutoff parameter (K max ). The values 9 (100-180 ba- 
sis functions per atom) and 10 (160-280 basis functions 
per atom) for R^^K max were used for the hybrid and 
semilocal calculations, respectively. The atomic sphere 
radius Rmt was chosen as 2.2 bohr for Rh and 2.1 bohr 
for Pd and Pt. 

The results of the calculations obtained with the LDA, 
PBE and YS-PBEO functionals are shown in Table Q] and 
Figs. [l][3] First, we start the discussion about the equi- 
librium geometry of the experimentally observed non- 
magnetic (NM) phase. The experimental values for ao 
(taken from Ref. |48| ) were extrapolated to T = K 
(using the volume expansion coefficient) and corrected 
for the zero-point oscillations according to the scheme 
of Alchagirov et al^ For Rh, we can see that the hybrid 
functional YS-PBEO gives a fairly accurate value of 3.799 
A (3.786 A for experiment), while LDA underestimates 
a by 0.03 A and PBE overestimates by 0.045 A. In the 
case of Pd, LDA, which leads to an underestimation of 
- 0.035 A, is more accurate than PBE and YS-PBEO 
which lead to rather large overestimations (> 0.06 A). 
For the 5d-TM Pt, LDA is more accurate (underestima- 
tion of only - 0.02 A) than PBE and YS-PBEO which 
yield too large equilibrium lattice constant. By compar- 
ing our values for ao to the PBE and HSE0G results from 
Ref. [l9| for Rh and Pd, we can see that the agreement is 
perfect (as expected) with PBE (less than 0.001 A of dif- 
ference), while our YS-PBEO lattice constants are about 
0.015 A larger than with HSE06. As explained in Ref. IH 
the differences between the YS-PBEO and HSE06 lattice 
constants are mainly due to the different schemes used to 
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FIG. 1: (Color online) (a) Total energy and (b) magnetic 
moment in the unit cell of the FM and NM phases of Rh ver- 
sus the lattice constant (expressed as the relative difference 
with respect to the experimental value) calculated from sev- 
eral functionals. For each functional, the zero of the energy 
was set to the lowest minimum. The vertical line indicates 
the experimental lattice constant. 

screen the semilocal exchange energy. In YS-PBEO, the 
simple scheme of Iikura et al— was used, while in HSE06 
the screening is done by considering the expression of the 
functional in terms of the exchange hole^i 

At this point it is useful to recall that the trends ob- 
served above with the semilocal functionals are actually 
rather general among the elemental TM (see, e.g., Ref. 
4). For the 3d series, PBE is the most accurate, while 
LDA strongly underestimates the lattice constants. For 
the id TM, GGA functionals which are softer than PBE 
(e.g., WC 33 or PBEsol 3 — ) are more accurate than LDA 
and PBE, while for the 5d, LDA is the most accurate 
method. It has also been observed (see, e.g., Ref. [l9l ) 
that for most solids (Li and Na are two exceptions), re- 
placing a fraction of semilocal exchange by HF exchange 
leads to a reduction of the equilibrium lattice constant, 
which is what has also been observed in the present work 



FIG. 2: (Color online) Same as Fig. [T]but for Pd. Note that 
in (a) and (b) the ranges for the x axis are different. 



with Rh, Pd, and Pt (i.e., the YS-PBEO lattice constants 
are smaller than the PBE ones). Therefore, the hybrid 
functionals based on PBE (HSE and YS-PBEO) correct 
the PBE overestimation for the lattice constantSi 13 i 19 
This also explains why the lattice constants calculated 
with the screened-exchange LDAii (100% of LDA ex- 
change replaced by HF exchange) are in most cases by 
far too small^ since on average LDA already under- 
estimates them. The results for the bulk modulus Bq 
(see Table HI show the usual trends: An underestimation 
(overestimation) of Bq is associated to an overestimation 
(underestimation) of clq. 

Spin-polarized calculations on Rh, Pd, and Pt were 
also done in attempts to stabilize a ferromagnetic (FM) 
solution. The results, whenever available, are shown in 
Table U and Figs. [T][3] In the case of Rh, it was not pos- 
sible to stabilize a solution with a non-zero magnetic mo- 
ment with LDA for the range of lattice constants that we 
considered (up to 15% larger than a cxpt ), while with PBE 
it is for values of a larger than a eX pt by about 13% (i.e., 
much larger than the equilibrium lattice constant) that 
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FIG. 3: (Color online) Same as Fig. [Qbut for Pt. Note that 
in (a) and (b) the ranges for the x axis are different. 



such a solution could be obtained. With YS-PBEO it was 
possible to stabilize a magnetic moment for lattice con- 
stants larger than the experimental one. Interestingly, 
the minimum of the FM total-energy curve is lower than 
the minimum of the NM curve by 0.59 mRy/atom [see 
Fig. [TJa) and Table |T]. This can be considered as a fail- 
ure of the screened hybrid functional which describes Rh 
as magnetic, and in addition the equilibrium lattice con- 
stant is more than 2% larger than in experiment. The 
value of the magnetic moment in the unit cell is around 
1.4 /j,b [see Fig. [ljb)] for lattice constants larger than 
a cxpt and drops suddenly to zero (NM solution) at this 
geometry. 

Turning now to Pd, which is known experimentally and 
theoretically to be on the verge of becoming FM,— we 
were able to stabilize a FM solution with all considered 
functionals. From Fig. HJb), we can see that the small- 
est lattice constant at which it was possible to obtain a 
FM solution is around —15%, —1.7% and 1% of a expt for 
YS-PBEO, PBE, and LDA, respectively, and therefore, at 
both the experimental geometry as well as at the corre- 



sponding theoretical geometry, the magnetic moment is 
non-zero for PBE and YS-PBEO. Note that the values of 
the magnetic moments obtained at larger values of a with 
the different functionals arc quite similar (0.40—0.45 /j,b)- 
As for Rh with YS-PBEO, the FM minimum is lower than 
the NM minimum for PBE and YS-PBEO functionals by 
0.15 and 3.96 mRy/atom, respectively. The latter value 
is one order of magnitude larger than in the case of Rh, 
and again this can be considered as a quite severe failure 
of the screened hybrid functional. 

For Pt [see Fig. [3Jb)], non-zero magnetic moments 
can be obtained for lattice constants larger than approx- 
imately -9%, 3%, and 7% of a oxpt for YS-PBEO, PBE, 
and LDA, respectively, and thus at the experimental or 
theoretical lattice constant only YS-PBEO gives a wrong 
magnetic ground state. Thus, from Fig. 02a) we can see 
again, that with YS-PBEO the FM state is more stable 
than the NM state and their minima differ substantially 
by 3.51 mRy/atom. 



III. SUMMARY 

From the results presented above, we can conclude 
that the screened hybrid functional YS-PBEO slightly im- 
proves over PBE for the equilibrium lattice constants of 
solids, which were too large with PBE (a rather general 
trend with PBE). But note that for the 4d-TM Pd and 
5d-TM Pt, LDA is more accurate than PBE and YS- 
PBEO. However, the main result of the present work 
is the inability of YS-PBEO to give the correct ground 
state for the three considered TM. Indeed, the FM state 
is found to be lower in energy than the experimentally 
observed NM state. It is well-known (see Ref. [l^ and ref- 
erences therein) that the use of the HF approximation for 
metals leads to severe problems like a vanishing density of 
states at the Fermi level or too large magnetic moments. 
(A remedy would be to use a compatible correlation func- 
tional like the random-phase approximation,— but such 
methods are very expensive and their self-consistent im- 
plementation not trivial at all.) 

These problems are in particular due to the long-range 
part of the HF exchange, which is not present in screened 
hybrid functionals. However, it was showni^ that the 
screened hybrid functionals still lead to an overestimation 
of the exchange splitting (and the corresponding mag- 
netic moment) in itinerant magnetic systems like Fe for 
which semilocal methods work quite well. 

While the use of the gradient correction was necessary 
to get correctly the body-centered cubic FM phase as the 
ground state in Fo 54 ' 55 (the ground state is face-centered 
cubic NM with LDA), our work has shown more evidence 
that using screened exact exchange as an additional in- 
gredient in functionals is not recommended for itinerant 
metals. This has important implications for the studies 
of correlated TM oxides on TM surfaces, which is rele- 
vant for many important catalytic active surfaces, where 
oxides on metals like Rh, Pd, or Pt play a crucial role. 
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While the hybrid functionals describe in such systems the Acknowledgments 
correlated oxides pretty well, the results suffer from the 
overestimation of magnetism for the underlying metals 56 
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